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Determining the heating mechanism (or mechanisms) that causes the outer atmo- 
sphere of the Sun, and many other stars, to reach temperatures orders of magnitude 
higher than their surface temperatures has long been a key problem. For decades 
the problem has been known as the coronal heating problem, but it is now clear that 
'coronal heating' cannot be treated or explained in isolation and that the heating 
of the whole solar atmosphere must be studied as a highly coupled system. The 
magnetic field of the star is known to play a key role, but, despite significant ad- 
vancements in solar telescopes, computing power and much greater understanding 
of theoretical mechanisms, the question of which mechanism or mechanisms are the 
dominant supplier of energy to the chromosphere and corona is still open. Follow- 
ing substantial recent progress, we consider the most likely contenders and discuss 
the key factors that have made, and still make, determining the actual (coronal) 
heating mechanism (or mechanisms) so difficult. 
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1. Introduction 

The coronal heating problem is one of the longest running solar physics puzzles 
and is still a highly controversial topic. Considerable progress has been made in 
modelling possible heating mechanisms, but discriminating amongst these to dis- 
cover if any one process dominates is extremely difficult to do. This is because the 
corona is not energetically isolated from the other regions of the atmosphere, such 
as the chromosphere, but instead the whole of the solar atmosphere forms a highly 
coupled system, with both energy and mass transferred in both directions between 
the chromosphere and corona through the transition region. The chromosphere has 
very different plasma properties to the corona and so explaining the heating of the 
solar atmosphere is now recognised to be considerably more complicated than had 
been appreciated during much of the twentieth century. 



The coronal heating problem first arose following the work of Grotrian ( 1939 ) 



(erroneously the author on the paper is given as W. Grotian) and Edlen ( 1942 ) who 
discovered that emission lines observed during a total solar eclipse in 1869 were 
not due to a new element called coronium, but were the result of highly ionised 
iron. This work demonstrated that the temperature of the corona is in excess of a 
million degrees Kelvin. In comparison, the temperature near the Sun's surface, the 
photosphere, is just 6,000 K. At the same time the density falls off by six orders of 
magnitude from the photosphere to the corona. 

The Sun is not unique in having such a hot outer atmosphere. Indeed, most stars 
in the cool half of the Hertzsprung-Russell (H-R) diagram have a corona and can 
be observed in X-rays. Their average temperatures range from 1-45 million degrees 
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Kelvin with the hottest associated with the most rapidly rotating stars. Using 



observations from the Einstein Observatory, Vaiana ( 1981 ) was the first to show 



that most stars, as opposed to just the odd exotic star, have soft X-ray emission 
and, hence, a corona. Linsky ( 1985 ) reviewed stars across the H-R diagram with 



emissions in ultraviolet and X-ray emission lines which are produced by plasma 
hotter than 10 4 K. He found that Dwarf stars of spectral type G-M and rapidly 
rotating subgiants and giants of spectral type F-K in spectroscopic binary systems 
are most likely to be solar-like in nature, i.e., they are likely to have "a turbulent 
magnetic field sufficiently strong to control the dynamics and energetics of their 
outer atmosphere". Furthermore, he found that T Tauri stars, other pre-Main- 
Sequence stars and single giants of spectral type F to early K are also probably 
solar-like. 





Figure 1. (a-e) SDO/AIA 211 A, 193 A, 094 A, 335 A and 304 A images of the solar corona 
and chromosphere showing the amazing complexity and vast range of scales of structures, 
(f) SDO/HMI magnetogram of the magnetic field in the photosphere. All images taken 
on 27th August 2011. 



X-ray, EUV and UV images of the Sun (Figure []J a-e)) show beautiful, finely 
detailed loops structures, as well as more diffuse emission which is likely to be 
formed from a mass of unresolved loops. The resolved loops cither connect regions of 
opposite-polarity magnetic field seen in the photosphere, or they are extended open 
structures that are anchored in a photospheric magnetic feature. As a general rule, 
the hotter the loops, the stronger the regions of magnetic field they are associated 



with (Fisher et al. 1998 1 



The structure of the magnetic field is highlighted in coronal images because 
heat is spread efficiently along field lines, but transport across them is greatly 
inhibited (thermal conduction in the corona is 13 orders of magnitude stronger 
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along field lines than across them (Priest 1982)). All of these factors signify that 



the magnetic field plays a major role in heating the solar corona. Moreover, there is 
strong evidence from other stars that X-ray coronae are associated with magnetic 
field since a power-law relation exists between X-ray luminosity and magnetic flux 
over many orders of magnitude, ranging from solar quiet-regions, through active 



regions, to dwarf and TTauri stars (Fisher et al. 1998[ Pevtsov et al.\ 2003 ). 



The loops observed, for example in Figure |l|a-e), are visible since they contain 
plasma of the appropriate density and temperature to emit at the given wavelength. 
The surrounding space contains loops that are either at a different temperature 
(hotter or cooler) or insufficiently dense to be observed. In equilibrium, and without 
any heating, a loop will be cool and rarefied, so any heating mechanism must explain 
both how the material fills the loop and heats it to coronal temperatures. 

In the corona, the combined radiative and conductive losses in active regions 

whereas in the 
_1 in active re- 



are 10 7 erg cm 2 s 1 and in the quiet-Sun are 3 x 10 5 erg cm 2 s 
chromosphere they are 2 x 10 7 erg cm~ 2 s _1 and 4 x 10 6 erg cm" 



gions and the quiet-Sun, respectively ( Withbroe & Noyes 1977). The source of this 
energy is convective churning of the plasma at and below the photosphere. These 
motions continuously move the footpoints of the ubiquitous magnetic fields that 
thread through the surface of the Sun with the timescale of these motions, relative 
to the end-to-end Alfven travel time, leading to one of two fundamental types of 
heating mechanism. Slow (long timescale) motions result in a quasi-static stressing 
of the field, whilst fast (short timescale) motions generate waves. The dissipation of 
magnetic stresses, which typically manifest themselves in the form of current sheets, 
leads to magnetic reconnection and is known as DC heating, whilst the dissipation 
of waves is referred to as AC heating. In both cases the actual dissipation will occur 
at kinetic scales, since the Lundquist number (the ratio of the Ohmic dissipation 
time to the Alfven crossing time) is very large (« 10 10 ~ 13 ). On the other-hand, 
the photospheric driving and the complex geometry/topology of the magnetic field 
are determined globally. Not surprisingly, this massive range of scales which need 
to be addressed by coronal heating models poses a problem for both theorists and 
observers. Current parallel computers have neither the memory nor the power to 
permit well resolved models on all scales and current imagining telescopes have 
neither the required spatial and temporal resolution nor the sensitivity to see any- 
thing at the kinetic scale. Additional difficulties arise in the form of the coupling 
between the dense interior and tenuous outer atmosphere which require very dif- 
ferent physics and, for open magnetic field regions, mechanisms must also account 
for the mass flux loss due to the solar wind. With all these difficulties it is not 
surprising that, to date, researchers have been unable to properly explain the solar 
heating problem. Nonetheless our understanding has advanced significantly. 

Since the identification of the hot solar atmosphere, there have been a huge num- 
ber of heating mechanisms proposed that involve reconnection or magnetic waves. 
A full review of the analytical and numerical studies carried out is far beyond the 
scope of this review, so we refer the interested reader to the reviews by Aschwanden| 
pOOll, iGoedbloed & Poedtsl pOOlh, iKlimchukl ((2006b, lAschwanden et al.\ 12007b, 



Goossens et al. (2011) and Reale (20101. Instead, we focus on those mechanisms 



currently most favoured. First, though, in Section 2, we consider the source of the 
corona's energy, before giving a brief selective discussion on the back grounds to 
wave heating (Section 3) and magnetic reconnection (Section 4). Nanoflare observa- 
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tions are reviewed in Section 5 and the heating of coronal loops, including numerical 
modelling and the recent observational results on spicules and Alfvenic waves, is 
the focus of Section 6. Finally, in Section 7, we present our conclusions. 

2. Poynting Flux 

The energy transported into or out of the solar atmosphere is provided by a Poynt- 
ing flux through the surface of the Sun, 



where E is the electric field, B the magnetic field, the magnetic permeability 
and S the solar surface. By assuming that motions are ideal and by writing the 
magnetic field and velocity in terms of components that are parallel or perpendicular 
to the solar surface (denoted by subscripts p or n, respectively) the Poynting flux 
per unit area through the solar surface can be written as 



where n is the normal to the solar surface directed radially outwards. The first 
term corresponds to the movement of existing magnetic footpoints by horizontal 
convective motions whereas the latter represents the emergence or submergence of 
magnetic flux. 

Typically the magnetic field is tilted back away from the direction of motion 
and so (B p • v p )B n /[i rts B p v p B n /fio provides an inflow of energy into the solar 
atmosphere (irrespective of the sign of B n ). If, instead the field was titled forward 
towards the direction of motion, then the Poynting flux would would instead remove 
energy from the solar atmosphere. Magnetic fields are thought to be held in intense 
kilo Gauss flux tubes and are moved about with a typical local speed of 1 km s _1 . 
Assuming that the magnetic field is inclined at an angle of 20 degrees, and that the 
filling factor is 1% in the quiet-Sun (giving a mean field strength of 10 G), then 
the Poynting flux per unit area is about 2.5 x 10 7 erg cm -2 s _1 . Active regions 
have much higher filling factors of around 50% (implying a mean field strength of 
500 G) which gives a Poynting flux per unit area of 1.2 x 10 9 erg cm~ 2 s _1 . In 
both cases simply considering the energy injected by moving the footpoints seems 
to give more than enough energy to heat both the chromosphere and corona, even 
if we assume that half the motions result in a loss rather than a gain in energy. Of 
course, though, energy may also be required to raise plasma from the chromosphere 
into the corona, accelerate particles and drive the solar wind in open field regions. 
With this in mind we also estimate the energy contribution from the emergence of 
flux. 

To obtain an inflow of energy to the solar atmosphere from the second term 
we require an upflow, and hence the emergence of flux. Newly emerging bipolar 
regions of flux become highly shredded as they rise through the top layers of the 
convection zone into the photosphere. They appear as a host of small-scale loops 
that are aligned so that after emergence they form a bipole with anything up to 
« 5 x 10 19 — 10 20 Mx of flux. Larger bipoles are found, but these are associated 




1 



(E x B) • n 



1 



(Bp-v p )B n + — (Bp-Bp)v. 
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with active regions. The Poynting flux associated with the emergence of a single 
kilo Gauss intense flux tube is 2 x 10 12 erg cm" 2 s" 1 , assuming an upflow of 2 km 
s -1 . However, the total em ergence rate of new flux is just 450 Mx cm -2 day _1 
(Thornton & Parnell 2011) suggesting a filling factor of 6 x 10~ 7 and so, in the 
quiet Sun, the Poynting flux due to the emergence of new magnetic flux is around 
10 6 erg cm~ 2 s _1 . This is a substantial fraction of the 4.3 x 10 6 erg cm~ 2 s" 1 that 
is required, but is a factor 25 less than that from the moving of footpoints. In active 
regions, the Poynting flux from moving flux is similarly dominant. 

The above provides a very simple order of magnitude argument, but its results 
are qualitatively in agreement with the results from the sophisticated MHD mod- 
elling of the upper convection zone to the lower corona using radiative transfer by 



Abbett & Fisher (2011 1. It is clear that the energy that heats the solar atmosphere 



comes from the convection zone and this energy manifests itself both as waves and, 
through the slow stressing, current sheets leading to reconnection. Which of these 
is most important for chromospheric and coronal heating is still very much an open 
question. 



3. Acoustic and MHD Waves 



Almost immediately after the identification in the early 1940's of highly ionised 
spectral lines, such as FelX and CaXIV, mechanisms based on the dissipation of 
acoustic waves were proposed (Biermann 1946 1948 Schwarzschildl 19481. How- 



ever, a heating theory based on waves must overcome several obstacles, namely: (i) 
waves must be generated in (or below) the solar surface layers, (ii) sufficient energy 
flux has to be transported as only a fraction of the wave energy will be transmitted 
into the corona and (iii) the waves have to dissipate (very) efficiently to convert 
the energy flux into heat in the right place, at the right time. The first of these 
aspects (i) is probably the easiest to address. Indeed, solar convection generates a 
mixture of upward propagating waves with an energy flux of several times 10 7 erg 
cm~ 2 s _1 (Narain & Ulmschneider 19961, which would be more than adequate to 



heat the solar corona (and accelerate the solar wind). However, (ii) is not so easy 
to satisfy since most of the (magneto) acoustic waves do not propagate into the 
corona due to strong reflection and refraction off the rapid density and tempera- 
ture increase in the Transition Region (in the traditional, plane-parallel view of the 
solar atmosphere). Additionally, mode coupling is expected to take place in regions 
where the sound and Alfven speed are approximately equal (see e.g. the numerical 



simulations by Rosenthal et al. (2002) and Bogdan et al. (2003)) reducing the flux 
reaching the corona still further. Alfven waves, on the other hand, are generally 
seen as very efficient carriers of energy over large distances ( Hollweg] 



1978 



1984 



Kudoh fc Shibata[|1999[ ). With regard to issue (iii), although [Wentzel| ( |1974| |1976[ ) 
pointed out that Alfven waves could heat the corona, their weak damping makes 
them problematic as a coronal heating mechanism. This led to the development 
of a variety of mechanisms which are likely to enhance the dissipation of Alfven 



waves such as resonant absorption (Ionson 1978 Rae & Roberts 1981 Goossens 



et al. 


1992 


Ofman et al. 


1998 


) and phase mixing ( 


Heyvaerts & Priest 


1983) 



Both resonant absorption and phase mixing are based on the property that in an 
inhomogeneous plasma, individual surfaces can oscillate with their own Alfven fre- 
quency In resonant absorption, on a specific magnetic surface, an incoming wave 
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motion can be in resonance with local oscillations, allowing the transfer of energy 
from large-scale motions to small lengthscales where dissipation can become ef- 
fective. Similarly, phase mixing describes how Alfven waves on neighbouring field 
lines (with different Alfven speeds due to a transverse inhomogeneity in the density 
and/or the magnetic field strength) gradually become out of phase, again generat- 
ing increasingly large transverse gradients until dissipative lengthscales are reached. 
Although mostly treated separately in the literature, examples of the close inter- 



play between resonant absorption and phase mixing can be found in e.g. Ruderman 



et al. ( 1997«|&" ). Apart from resonant absorption and phase mixing, other mecha- 
nisms such as (non-linear) mode conversion (to modes which might more readily 
dissipate) and ion-cyclotron resonance have also been proposed (see for example, re- 
Aschwanden| ( [2004| ; |Goedbloed fc Poedtslp^jKlnrTcliTu^ 



views by 



& Ballai 



( |2007 ) and references therein). 



Although many wave-based heating theories were developed, especially in the 
1980 'sand early 1990's, (See|Aschwanden||2004|[Goedbloed fc Poedts[[2004l|Goossen|s 



et al. 2011, for more detailed reviews) reconnection-based mechanisms were gener- 



ally favoured to address the coronal heating problem. Porter et al. (1994) point out 
that observed non-thermal line broadenings imply sufficient wave fluxes to account 
for the energy requirements of both quiet Sun and active regions, but that this does 
not necessarily mean that these fluxes are readily converted into heat. Indeed, even 
if the line broadenings are actually produced by waves, the efficiency of many of the 
theoretically studied wave-based heating mechanism relies on the presence of waves 
of the right type and the right frequency, which, at the time, could not be verified 
since direct observations of waves in the solar atmosphere were largely absent. This 
situation has changed with the advent of space-based EUV imagers (SOHO/EIT 
and TRACE), which now have sufficiently high spatial and temporal resolution to 
resolve oscillatory motions. A large variety of waves and oscillations have since been 
observed in the solar atmosphere. However, the observed energy budget contained 
within the observed waves was generally several orders of magnitude too small to 
account for the heating of coronal loops. Additionally, given the rapid damping of 
many of the observed oscillations, it is not at all clear that the timescales involved 
even allow for heating to occur during the actual oscillations. Hence, wave-based 
coronal heating mechanisms initially gained little ground. The discovery of several 
types of MHD waves in the solar coronal did however led to the rapid development 
of coronal seismology (see the review by De Moortel & Nakariakov 2012 in this 
issue). 

More recently observations of flows and Alfven waves have led to a renewed 
interest in wave heating mechanisms (see Section 6c and 6d), both observationally 
and theoretically. For example, Van Doorsselaere et al. ( 2007 1 look at coronal loop 
heating profiles and find that, based on observational evidence, a resistive wave 
heating mechanism is more likely than a viscous one. Similarly, Taroyan et al. ( 2007 I 
show that power spectra of Doppler shift oscillations might allow a distinction 
between models based on footpoint heating and uniform heating. A recent study 
by Antolin & Shibata (2010) investigates heating based on torsional Alfven waves 
and suggests that, within the constraints of their 1.5D model, this mechanism is 
only viable in a rather narrow regime, requiring long and thick loops 



(sec 



Jess et al. 



(2009) for a possible observational identification of torsional Alfven waves in the 
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solar chromosphere). We refer the interested reader to Taroyan & Erdelyi (2009) 
for a review of heating diagnostics with MHD waves. 



4. Magnetic Reconnection 

One consequence of new flux (> 10 19 Mx) emerging into the photosphere is the 
formation of new coronal loops connecting the new flux features to pre-existing ones. 
The only way these new loops could have formed is through magnetic reconnection 



(see article by Pontin in this volume or books by Priest & Forbes (2000) and 
Biskamp (2000) for more details) allowing a global restructuring of the magnetic 
field. The removal of flux from the photosphere is also associated with magnetic 
reconnection since previously unconnected pairs of magnetic features must connect 
before they can cancel. In both cases, the reconnection driven by emergence or 
cancellation of flux may be associated with not only the formation and heating 
of new loops, but also the bulk acceleration of plasma, particle acceleration and 



the lifting of mass into the chromosphere/corona (e.g. Archontis et al. 2005 



Rekowski & Hood 2008 Martinez- Sykora et al. 2009 ). Depending on the location 



von 



and rate of reconnection and also the amount of plasma lifted, such processes may 
result in loops that can range from hot coronal loops to cool dense chromospheric 
loops. 

Emerging and/or cancelling flux models have been used to explain X-ray bright 



points ( 


Heyvaerts et al. 


1977 


Priest et al. 


1994 


Parnell et al. 


1994 


Longcope 


1998 


) and X-ray jets ( 


Shibata et al. 


1992 a\b 


Yokoyama & Shibata 


1995 


Moreno- 



Insertis et al. 2008), although, the basic principles of these models also work on 



larger scales in active 


regions (e.g. 


Masuda et al] 


1994 


Archontis et al. 


Martmez-Sykora et al. 




2009| Parnell et al. 




2010) 


indicating that many 



2005 



observed coronal loops, large and small, may well result from reconnection. 

The processes of emergence and cancellation, however, are not essential for re- 
connection to occur. The mere movement of magnetic footpoints can also lead to re- 



connection ( 


Longcope 


1998 


Galsgaard et al. 2000 Aulanier et al. 


2006 


De Moor- 


tel & Galsgaard 


200G 


Haynes et al. 


2007 


Aulanier et al. 


2007 


)■ 


Close et al. 


(2004) 



used a magnetic feature tracking algorithm to identify and track magnetic features 
observed in quiet-Sun MDI high-resolution magnetograms. They then extrapolated 
the magnetic field into the solar atmosphere and quantified the changes in mag- 
netic connections between features as the photospheric magnetic field evolved. Their 
work suggested that all the flux in the corona is recycled (all the connections are 
broken and reconnected) in around 1.4 hrs, a rate at least ten times faster than 



the recycling time of the photospheric flux as observed by MDI (Hagenaar 2001 ). 
Furthermore, since the flux in the photosphere is contained predominantly within 
small-scale magnetic features, the complexity in the structure of the magnetic field 
falls off rapidly, with more than 50% of the flux closing down below 2.5 Mm and 
just 10% reaching above 25 Mm (Close et al.\ 2003). These results indicate that 



the majority of the magnetic interactions (reconnections) between loops will occur 
below these heights, quite probably in the chromosphere. 



The models of ( Close et al. 2003 2004 ) are oversimplified and more sophisti 



cated models are required to quantify these effects more accurately. For instance, 
in the solar atmosphere reconnection occurs more readily in the cool dense chromo- 
spheric region than in the rarefied hot corona since diffusivity is inversely propor- 
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tional to temperature. Furthermore, reconnection will convert magnetic energy into 
both thermal and kinetic energy (generating waves and shocks) , as well as particle 
acceleration. In particular, there will be direct heating from ohmic dissipation of 
the plasma in field lines threading the reconnection site, as well as indirect heating 
of the surrounding plasma due to the viscous damping of waves, outflow jets and 
shocks created by the reconnection. In some reconnection models, e.g., |Petschek| 
(1964) it is believed that the majority of the heating comes from the shocks. Al- 
though, recent numerical experiments suggests ohmic heating dominates, in both 



the chromosphere and corona, for some reconnection scenarios (Fuentes Fernandez 



et al. 


2012 


Fuentes Fernandez & Parnell 


2012) 



Thus it is still an open question as 
to whether viscous or Ohmic heating mechanisms dominate and what the location 
of the energy deposition is relative to the reconnection site. However, it is clear that 
the precise nature will depend on the particular reconnection regime operating and 
the properties surrounding plasma. 

Gudiksen fc Nordlund] ( |2005a|& 1 were the first to run 3D MHD models solving an 



energy equation that included field aligned thermal conduction, as well as radiative 
losses. Their results qualitatively match coronal observations in many ways, but in 
such complex numerical experiments it was not possible, at the time, to identify 
exactly which physical processes were responsible for the heating. For instance, 
some key questions which need to be addressed are, in their model is reconnection 
the key heating mechanism and if so is the reconnection occurring principally in the 
corona or the chromosphere? Also, the cadence of their experiments was such that 
analysing the role of waves was not possible. Yet more sophisticated modelling (e.g., 



Martrnez-Sykora et al. 2009 Hansteen et al. 2010), including additional physics 
such as optically thick radiation in the chromosphere, has also been performed. 
These models find episodic heating events in the chromosphere and it is claimed that 
these events are due to 'stressing of the magnetic field' by photospheric motions, 
as well as flux emergence, but a detailed analysis on the true nature of these events 
has yet to be undertaken. 



5. Nanoflares - Observations 



One important heating mechanism, the basics of which were proposed almost forty 
years ago, is the dissipation of a multitude of small-scale current sheets (Parker 



1972 Levine 1974 1. In particular, Parker (1972) hypothesised that slow convective 



motions of magnetic foot points could lead to a braiding of the magnetic field, 
creating small current sheets which would be continually dissipated and reformed, 
providing direct heating via Ohmic dissipation. This concept may be thought of 
as a braiding between distinct magnetic loops or a self braiding of the magnetic 
field within a single loop. In either case, magnetic loops need not necessarily wrap 
around each other, but simple shearing motions suffice to produce current sheets. 



Following the discovery by Lin et al. ( 1984 ) from rocket-borne observations 
of many small hard X-ray spikes with energies between 10 24 — 10 27 ergs, Parker 



(1988) speculated that a sea of small flares ('nanoflares'), with energies nine orders 
of magnitude less than large flares (i.e. 10 24 ergs), could heat the corona if they 
occurred in sufficient numbers. This is known as Parker's braiding or nanoflare 
theory. 
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In reality a wide range of events (sizes and energies) have been observed in 
the solar atmosphere following power-law distributions of frequency with varying 
indices. Hudson ( 1991 ) considered the implications of a power-law distribution of 



energies. His important paper pointed out that the index of the power-law distri- 
bution of energies is critical in establishing whether small-scale or large-scale flares 
provide more heat. If the powerlaw index a > — 2, large flares dominate the heat- 
ing of the corona, whereas if a < —2, small flares dominate. Note, though, even if 
the slope is steep enough, nanoflares can only heat the corona provided there are 
enough of them to produce sufficient energy. 

Using HXRBS/SMM, Yohkoh/SXT, SOHO/EIT and TRACE t he distribution 
of flares, microflares and/or nanoflares were determined by m any (jCrosby et al. 



1993 



et al. 



Shimizu| p95l |Krucker fc Benz[ [19981 |Parnell fc Jupp[ |2000[ |Aschwandeh 



2000] [Benz & Krucker[ |2002| |Aschwanden & Parnell[ |2002| . Power-law indices 



for the distribution of event energies were reported ranging from —2.7 < a < —1.5 
leaving the validity of nanoflares as a heating mechanism completely open. Similarly, 
power-law indices in the range —2.5 < a < —1.85 have been reported for stellar 



flares observed on a range of late-type cool stars (spectral type F-M) ( Audard et al. 



2000 Gudel et al. 2003 Arzner et al.\ 2007). 



Observational determinations of a come from the analysis of images taken by 
different instruments (e.g., SOHO/EIT, TRACE), with different detection routines 
(e.g., local peaks in emission or intensity in particular lines in a specific order of 
some specified size in area and peak), different energies (e.g., thermal energy or 
radiative/conductive losses) determined and different assumptions made. Parnell 



(2004) examined the effects of making different, but reasonable, assumptions to 
detect flares, determine their energies and calculate the resulting slope of the power 
law. The approach involved using just two data sets, one of the quiet Sun and one 
of an active Sun taken by TRACE. 

In the end Parnell (2004) determined 1200 different estimates of the power-law 
index by considering different sets of the possible parameters as used by others 
in their analysis. For both data sets (quiet-Sun and active-region) a lay in the 
range —2.5 to —1.6 with a mean of —1.9. The most reasonable conclusion that one 
can draw from these results is that it is not possible to accurately determine the 
distribution of individual nanoflarc/flare energies using the method of direct detec- 
tion. Instead, alternative methods that attempt to detect ensembles of nanoflares 
have been sought, (e.g. Sakamoto et al. 2008 Terzo et al. 2011) providing sup- 



port, but not conclusive observational proof, that nanoflares heat the corona. Thus, 
nanoflare heating is currently still regarded as one of the most promising coronal 
heating mechanisms for both coronal loops and the background quiet-Sun. 



6. Heating of Coronal Loops 

Figure [l] clearly shows that the solar corona is made up of many loops which range 
from just a few Mm in length to hundreds of Mm. Even the diffuse emission seen 
in many of the images in Figure [T] is highly likely to be made up of unresolved 
mini loops. The footpoints of coronal loops are not made up of single monolithic 
magnetic features of one polarity, but are made up of many tiny intense flux tubes, 
or strands, which are continuously being sheared, rotated and compressed as a 
result of the underlying granular and supergranular motions. Modelling has shown 
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that any one of these motions can create tiny current sheets within a coronal loop, 
i.e., as in Parker's braiding/nanoflare (Parker 1988), the kink instability model of 
Hood et al. (2009) or the tectonics idea of Priest et al. (2002). The idea is that, if 
there are enough current sheets distributed throughout the loop cross-section, and 
these are dissipated and regenerated at an appropriate rate, then enough energy 
may be provided to maintain the loop at more than a million K for several hours, as 
observed. However, it is also possible that small-scale heating events may result from 



a range of different wave dissipation mechanisms such as phase mixing (Heyvaerts 



& Priest, 1983) and resonant absorption (Ionson 1978 Hollweg & Yang 



1988). 



(a) 3D MHD Loop Modelling 



Models in which simple magnetic fields have been shown to produce multiple 
current sheets through random driving motions have been conducted by many (e.g., 



Mikic et al. 1989 Galsgaard fc Nordlund 1996 Rappazzo et al. 2007 2010). For 



instance, Galsgaard & Nordlund ( 1996 ) conducted 3D MHD numerical experiments 
in which they shear an initially uniform vertical magnetic field first in one direction 
then at right angles. The driving pattern was repeated continuously with random 
driving periods in each direction forming a quasi-steady state, involving a chaotic 
distribution of multiple small-scale current sheets throughout the domain where 
the current sheets continuously dissipate and reform elsewhere. Each small-scale 
reconnection event may be thought of as a single nanoflare which locally heats 
a few strands of the coronal loop. However, by averaging over the depth of the 
domain and over a short time period, the loop appeared to be approximately uni- 
formly heated (Figure If the same braiding experiment is carried out, but 
this time with a stratified loop, thermal conduction and optically thin radiation, 
then within the chromospheric regions of the loop considerably more current sheets 
are both formed and dissipated in comparison to the higher temperature coronal 
regions (Figure [2jb)) (Gals gaard] 2002). This imbalance in the distribution of heat- 
ing events is due to the significant (a factor of 10) increase in Alfven speed from the 
photosphere to the corona. The energy density from propagating perturbations is 
therefore much higher, per unit length, below the transition region than above. Fur- 
thermore, the transition region acts as a reflector preventing much of the Poynting 
flux injected by the boundary motions from reaching the corona. 

Instead of a complex driving pattern on a simple magnetic field creating multiple 
forced reconnection events a multitude of small-scale reconnection events can also 



be created from a complex magnetic field which underg 


oes an instability (e.g., 


Hood 


et al.\ 


2009 


Wilmot-Smith et al. 


2010 


Pontin et al. 




2011 


)• 


Wilmot-Smith et al. 



(20101) and Pontin et al. (20111 considered the resistive relaxation of a coronal 



loop consisting of a complex braided magnetic field in which an instability leads to 
spontaneous reconnection producing a multitude of tiny nanoflares (Figure [2|c-e)). 
They conducted 3D MHD experiments in which this initially braided magnetic field 
(containing just 6 crossings), formed by slow footpoint motions, is allowed to relax 
via resistive MHD. Initially, the magnetic field has large-scale diffuse currents, but, 
as the field evolves, intense current layers are naturally generated until a resistive 
instability spontaneously triggers reconnection and the start of a cascade in energy 
release to small scales. Most importantly they show that in a magnetic loop with 
what appears to be just 6 reconnection sites (braidings), many hundreds of little 
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Figure 2. The deposition of energy vs position along a curved loop due to the braiding 
of a uniform magnetic field in a loop with a (a) uniform atmosphere and (b) a stratified 
atmosphere (from Galsgaard ( 2002 1). (c) Sample magnetic field lines showing the initial 
force-free field in the braiding experiment by Pontin et al. (2011 1, (d) and (e) show isosur- 
faces at the start and part way through once reconnection is underway showing how the 
large current structures have fragmented into a multitude of small-scale sheets providing 



widespread heating throughout the loop (from Pontin et al. ( 2011 1) 



energy release events occur, distributed throughout the coronal loop (Figure [2jc)). 
Moreover, there is much more reconnection than one would first predict. This is 
because the magnetic field is not clever: it does not follow the route to the quickest 
and simplest untangling of itself, but rather reconnects wherever the conditions for 



the onset of reconnection are met (Parnell & Haynes 2010 Parnell et al. 2011). 



This results in a small reduction in the magnetic stresses in the field, but not a 
complete untangling, and a continuous chain of reconnection events must occur 



before the field can completely untangle itself. Indeed, Pontin et al. (20111 find 



that the flux in their braided flux tubes can reconnect up to 4.4 times, confirming 
that the flux is not untangled via the quickest route. Such a process is known as 



multiple (or recursive) reconnection and was first found by Parnell et al. (2008) 



With so many small-scale energy releases throughout the coronal loop, the entire 
loop can easily be heated and not just a few strands from within it. 



Parnell & Haynes (20101 showed that multiple reconnection has a series of im- 



portant consequences for coronal heating. Not only does it result in a multitude 
of small-scale reconnection events, and therefore widespread heating throughout 
the coronal loop, but considerable energy maybe released over a prolonged period. 



Parnell & Haynes (2010) considered an experiment where two opposite-polarity 



magnetic sources, in a background uniform overlying field, are driven past each 
other at a constant slow speed of 5% of the Alfven crossing time. They ran a series 
of numerical 3D MHD experiments which were exactly the same except for the 
value of the resistivity. They found that as the resistivity decreased and became 
more realistic the magnetic reconnection between the opposite-polarity sources was 
more complex in nature and thus took longer. This resulted in an increase in the 
horizontal component of the magnetic field at the base and thus an increase in the 
Poynting flux injected into the system. With more energy in the system, more en- 
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Figure 3. (a) Graph depicting the temporal evolution of the temperature (dot-dashed), 
light curve (solid) and density (dashed) in a loop(strand) due to pulsed heating from 
a single coronal heating event, (b) A temperature vs density graph of a loop during the 
same event (Phases I- IV). A long-duration heating pulse may lead to an evolution through 
phases I, II* and IV*. The dotted line indicates the equilibrium state in which conductive 



and radiative cooling balance. Graphs based on Figures 14 and 15 from Reale (2010). 



ergy had to be released to get back to the same relaxed state in each case. In their 



braiding experiment Pontin et al. (2011) note similar findings. Thus, the nature, 



location and timing of magnetic energy release events are important to understand, 
not only because of their feedback on the system which may affect the amount 
of Poynting flux subsequently injected, but also because the plasma/observational 
consequences will depend greatly on these factors too. 



(b) ID Loop Modelling 



Numerical 3D MHD experiments, such as those described above, are very useful 
in that they can give a good idea as to the distribution, frequency and size of 
energy release events and can provide an indication of the mcchanism(s) heating 
the corona. However, due to the large global scales they cover they rarely include 
all the appropriate physics. In particular, a very simple energy equation is often 
used neglecting thermal conduction, radiation, radiative transfer and ionization 
non-equilibrium effects. Thus the plasma response seen in MHD models is not 
correct and so comparing with observations is extremely difficult. However, the 
very limited energy transport across fieldlines means that ID hydro models which 
can incorporate many of these important physical processes can be used to model 
the plasma response in a single strand, and a whole coronal loop, by combining 
the results from many single strand runs (see Klimchuk 2006 Reale 2010 for a 
review) . 

It is highly unlikely that cool dense loops thread the corona waiting to be heated. 
Instead any coronal heating mechanism must explain both how the plasma is trans- 
ported into the corona and how it is heated. One crucial result of ID hydrodynamic 
models has been to create a clear picture of the behaviour of the plasma within 
a strand during an impulsive coronal heating event (note that, although often re- 
ferred to as a 'nanofiare', this event may be either AC or DC in origin) from the 
starting point of an initially cool, low density, equilibrium strand, i.e., one that 



will not be emitting in EUV or X-rays ( Cargill & Klimchuk 2004 Klimchuk 2006 



Reale 2007). An impulsive coronal event will cause the strand to evolve through 



four phases (Figure [3]) . During the first phase (I) , heat from the event is conducted 
along the field rapidly raising the temperature of the whole strand up to a max- 
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imum. The density and emission show very little change during this very short 
phase creating a strand that is underdense relative to static equilibrium theory. 
A second longer phase (II), during which the impulsive event continues to heat, 
follows in which the temperature remains fairly constant, but the chromosphere is 
strongly heated causing it to expand leading to the evaporation of plasma which 
rises and fills the strand. During this so-called evaporation phase, the emission rises 
rapidly and the strand becomes denser, but it still remains underdense. The third 
phase (III) starts when the impulsive event stops and conduction starts to cool the 
plasma. The emission first peaks, then starts to decrease, but the density continues 
to increase. As the strand cools and becomes more dense, radiative cooling effects 
becomes more important until, at the end of this phase, radiative and conductive 
cooling finally balance. In phase IV, the loop appears overdense, radiative cooling 
dominates and the density starts to decrease, slowly at first, until the temperature 
decreases to a point at which catastrophic cooling occurs and the density plum- 
mets until the temperature, density and emission measure return to their pre-event 
values. 

Clearly, this scenario has a number of key observables that can be tested. 
Namely, at the start, extremely hot plasma, which, it has been calculated, must 
be in excess of 4 MK, should be present (Cargill & Klimchuk 19971. However, 



since the loop plasma density is low during this phase, any emission will be very 
weak and so such hot plasma will be difficult to detect. Schmelz et al. (2009) have 
observed hot plasma outside solar flares which is widespread throughout active re- 
gions, but weak. A recent paper by Reale et al. (2011) claims to have now found 



thin high temperature strands of around 10 MK within coronal active regions us- 



ing SDO/AIA. However, such hot emission is likely to be rare since Bradshaw & 



Klimchuk (2011) argue that noncquilibrium ionization effects will suppress the hot 



emission lines and in some cases render them invisible. 

Secondly, evaporation upflows should be present. There appears to be a straight- 
forward split with hot chromospheric UV emission lines and coronal EUV lines, in 
the quiet Sun, with temperatures of formation above 0.5 MK having blue Doppler 
shifts (e.g., Ne VIII, Mg X, Fe XII), whilst the cooler (< 0.5 MK) emission lines 
(e.g., C II, C IV, Si IV, S VI, N V, O IV, O V, O VI) are red shifted ([Peter fc Judge 



1999). In active regions the split is normally higher at around 1 MK. These results 



may seem consistent with a hot evaporative upflow and a cooler downflow, how- 
ever, there are other interpretations. In particular, the hot upflows may very well 
indicate that the actual heating mechanism is not located in the corona, but that 
plasma in the chromosphere and transition region is heated in situ, resulting in its 
expansion, creating a flow which fills the loop with hot plasma. There is currently 
a hot debate over which interpretation is correct, especially since a direct signature 
of evaporation such as the impact of accelerated particles on the chromosphere is 
likely to be weak and, hence, very difficult to identify. Most likely though both 
mechanisms are in operation, but which is most dominant is still an open question. 

Finally, since the density behaviour lags that of the temperature during the 
heating, evaporation and conductive cooling phases (I, II & III) the strand appears 
underdense relative to static equilibrium predictions. Observations suggest that 
hot loops (T > 2 MK) are underdense (Porter & Klimchuk 1995), and so they are 
consistent with loops made up of multiple impulsively heated strands where the 
majority of strands are in Phases I, II and III. However, so-called warm EUV loops 
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(T w 1 — 2 MK) are observed to be overdense relative to equilibrium models and are 
found to have a narrow temperature range, appearing practically isothermal along 



their length (Lenz et al. 1999 Winebarger et al. 2003 ). More recently, observations 



from EIS indicate that the plasma in these loops is not actually isothermal, but 
rather multithermal with a narrow temperature range (10 6 ± 3 x 10 5 K) and a 
volumetric filling factors of just 10% (Warren et al. 2008[ ). Thus it is thought that 
these warm loops are highly filamentary and cooling from a hotter state (Warren 



et al. 2002 2003 ) . A possible explanation is that the loop is made up of multiple 
strands each of which is impulsively heating by a single, short-lived (in comparison 
to the characteristic cooling time) event. Hence, the plasma in the strand evolves as 
it is first heated and then cools, such that overall the loop appears to be overdense 
(in phase IV) for an extended period (KhmchukJ 2006). 

For a heating event that is longer than a conductive cooling time the plasma in a 
strand may follow the phases I, II* and IV* (Figure^) instead. In fact strands that 
are heated by either a long slow heating event or a series of repeated short impulsive 
events could evaporate sufficient density such that the strand reaches a hydrostatic 
equilibrium and may remain in equilibrium (between phases II* and IV*) if the 
heating balances the conductive and radiative cooling. When the heating event 
ceases, radiation is likely to dominate and the temperature and density decrease 
(Phase IV*) ( |Klimchuk et aL[|2006||Reale[|2007| ). Some loops have a wide range of 
densities and lengths and, in general, follow a series of scaling laws derived from 



hydrodynamic equilibrium modelling (such as Rosner et al. 1978). They appear 



to be consistent with steady heating, i.e., a deposition of heating at a rate that 



balances the conductive and radiative cooling of the loop (e.g. Porter & Klimchuk 



1995 Kano & Tsuneta 1995 Warren et al. 2010 1. Indeed, models of entire active 



regions or the whole global corona have been made by applying these scaling laws 



to potential and non-linear force-free magnetic field extrapolations ( Schrijver et al. 


2004 


Warren & Winebarger 


2006 Winebarger et al. 


2008 Lundquist et al. 


2008) 



The coronally heated, multistrand scenarios discussed above fit quite a few, but 
not all observations. In particular, their main problem is that they do not explain 
how the chromosphere itself is heated, although, if the heating events are due to 



braiding, then, as the 3D braiding experiments of Galsgaard ( 2002 ) seem to suggest, 
this should not be a problem as an even greater number of heating events should 
occur in the chromosphere than in the corona. Furthermore, as discussed in the 
introduction, waves excited by footpoint motions are likely to be abundant in the 
chromosphere. Nevertheless, a recent series of observations and numerical experi- 
ments now suggest that the corona may actually be heated by strong upflows that 
are driven from the photosphere/chromosphere, as discussed in the next section. 



(c) Spicules & Propagating Disturbances 
Following the launch of SOHO and TRACE, several authors such as Bcrgh mans| 



fc Clette| ( |1999| ) and |De Moortel et aL| ( |2000[ ), reported the observat ion of propa- 



gating disturbances in large, quiescent (fan) loops at the edges of active regions, 
which were interpreted as slow magnetoacoustic waves propagating along coronal 
loops. However, the energy budget contained within these observed propagating 
disturbances was generally several orders of magnitude too small to account for the 
heating of coronal loops. Even though the coronal energy budget of the propagat- 
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Figure 4. (left) An SDO/AIA 171 A snapshot of the active region analysed by Mcintosh 
(2011). Inset is an unsharp-masked image of the region outlined in the red box, 



et al. 



showing the transverse loop(thread) displacements, (right) Schematic representation of the 
footpoint driven Alfvenic perturbations in open and closed field geometries (representing 
coronal holes and quiet-Sun/active-region loops, respectively). The red arrows represent 
the oscillatory displacements observed by |McIntosh et al.\ ( |2011[ ) in coronal loops. 



ing slow waves was thought to be very small, Fontenla et al. (1993) and Jefferies 



et al. (2006) point out that in the lower atmosphere, the energy flux contained in 



the low- frequency waves (< 5 mHz) is much higher than the flux in their high- 
frequency counterparts. These high-frequency waves had long been considered as 
prime candidates to address the heating of the chromosphere, an idea which was 



rejected by Fossum & Carlsson (20051. However, their low- frequency counterparts 



seem to contain much more energy than previously thought. 

The potential of the low-frequency propagating disturbances has very much 
come to the forefront again with recent observational results from Hinode and 
SDO/AIA as they have been linked with the supply of mass and energy to the solar 



corona and solar wind. De Pontieu et al. (2009) and De Pontieu et al. (2011 ) suggest 



that a significant fraction of the coronal mass is supplied by chromospheric jets 
("Type II spicules"). These authors analyse Hinode and SDO data, in which they 
find faint, but ubiquitous, quasi-periodic disturbances, interpreted by these authors 



as flows rather than slow magnetoacoustic waves (see De Pontieu & Mcintosh 



2010 Verwichte et al. 2010 for more details on this debate), and estimate that 



the mass contained within these spicules could supply a substantial part of the 
coronal hot plasma. Recently, Mcintosh et aL| (2012) reported on observations of 
slow, countcrstreaming downflows (~10 km/s) in cooler spectroscopic lines, which 
the authors suggest could be the return-flow of coronal material (i.e. the end of the 
mass cycle). Although this mechanism is under considerable debate, it certainly 
highlights that the corona cannot be treated in isolation. 



(d) Transverse, Alfvenic Waves 



As well as quasi-periodic propagating disturbances (flows/waves) propagating 
along coronal loops, recent observations have revealed transverse, propagating os- 



cillations along a variety of structures such as spicules ( De Pontieu et al. 2007 



He et al. 2009a|& ), X-ray jets (Cirtain et al. 2007), prominence threads (Okamoto 



et al. 20071 and coronal loops (Tomczyk et al. 2007 Tomczyk fc Mcintosh 2009 



Mcintosh et al. 2011 ). An example of the active region analysed by Mcintosh et al. 
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(2011) is shown in Fig. |4j who report on small-amplitude transverse displacements 
in coronal holes, quiet-Sun and active-region loops. An example of the loop(thread) 
displacements can be seen in the unsharp masked image inset in Fig. [4] (left). As 
these footpoint-driven, oscillatory displacements are interpreted as 'Alfvenic' (see 
schematic presentati on in Fig.[4|, they are of particular interest to the coronal heat- 
ing question. Indeed, De Pontieu et al. ( 2007[ ) estimate that the energy flux carried 
by chromospheric Alfvenic waves (observed by Hinode/SOT) is sufficient to accel- 
erate the solar wind and possibly heat the quiet-Sun corona. Although |Tomczyk| 



et al. (2007) report that the energy budget estimated from the COMP observations 



of propagating disturbances in off-limb coronal loops is too small to account for 
coronal loop heating requirements, Mcintosh et al. (20111 not only demonstrate 



that small-amplitude, transverse, Alfvenic perturbations are ubiquitous in the so- 
lar corona but that the energy contained in these oscillations appears sufficient 
to address the requirements of the quiet Sun and coronal holes (~ 10 5 — 2 x 10 5 
erg cm~ 2 s _1 ). In active regions, the reported energy (~ 10 5 erg cm -2 s _1 ) is too 



small to acco unt for the heating of coronal loops (which according to Withbroe & 
Noyes ( 19771 would require about 2 x 10 6 erg cm~ 2 s _1 ). However, in active regions, 



the superposition of a large number of randomly directed oscillating loop(strands) 
along the line of sight implies that the estimate reported by Mcintosh et al (2011) 
is likely to be a lower limit. Similarly, a large number of loops, with transverse 
displacement in random directions is likely to be present off-limb along the line 
of sight, which could explain why the energy budget estimated by Tomczyk et al 



(2007) seemed very small. Using 3D MHD numerical simulations, De Moortel & 



Pascoe (2012) show that the energy budgets estimated from observed, line-of-sight 
integrated Doppler velocities (such as those observed by Tomczyk et al. ( 2007| )) will 
be at least an order of magnitude smaller than the energy present in the 3D coronal 
volume. 

Although there is currently still considerable discussion about the interpretation 
of the propagating disturbances and the transverse displacements, what is clear is 
that the role of waves in the coronal heating problem has to be reassessed in light of 
the recent observations by SDO and Hinode. It now appears that Alfvenic waves are 
ubiquitous in the solar atmosphere, potentially containing a substantial amount of 
energy, and that chromospheric jets form an integral part of the mass supply cycle 
for the corona and solar wind. However, this does not necessarily imply a solution 
to the coronal heating problem. Although observations have unambiguously shown 
that waves are present, theoretical modelling now needs to address how these lead 
to heating, in the right locations, at the right time. Indeed, one has to keep in mind 
that although there appears to be a considerable amount of energy in the observed 
Alfvenic waves (see Mcintosh et al. 2011 ), this does not translate automatically into 
heating of the plasma. In particular, several authors have recently interpreted the 
observed damping of, for example, the transverse oscillations observed by COMP 
(Tomczyk et al. 2007 Tomczyk & Mcintosh 2009) as mode coupling, whereby 



the amplitude of the footpoint driven, observed, transverse oscillations decreases 



due to mode coupling to an azimuthal Alfven wave in the tube boundary (see De 



Moortel & Nakariakov 2012, for a detailed discussion). Potentially, the dissipation 
of the azimuthal Alfven waves, enhanced by phase mixing in the inhomogeneous 
tube boundaries could lead to local plasma heating. However, the viability of this 
remains to be demonstrated. 
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7. Discussion Sz Conclusions 

For decades after the discovery of the very high temperatures in the solar atmo- 
sphere, solving the coronal heating problem was mostly framed as an 'AC or DC 
question', or 'waves versus reconnection'. Within the last few years, observational 
results have shown that the coronal heating problem has to be thought of in a dif- 
ferent way. The highly dynamical nature of the solar atmosphere will clearly lead to 
the occurrence of both types of processes and the true question is to determine the 
relative contribution of the different mechanisms in different coronal structures. 
Wave-based heating mechanisms very much seem to have come full circle, from 
being the first mechanism suggested to provide the necessary energy input into 
the solar atmosphere, to being mostly discarded and recently back to the forefront 
again, driven by the spectacular observations by Hinode and SDO. Additionally, 
the discovery of the abundance of Type II spicules demands that the impact of the 
equally dynamic chromosphere has to be assessed much more carefully and that 
it is essential that the 'coronal' heating question must be studied as part of the 
highly coupled solar atmosphere. Indeed, it is well known (but often ignored) that 
the energy requirements of the chromosphere far exceed those of the corona. 

Heating by multiple small-scale reconnection events has long been considered a 
viable mechanism and much support has been lent to this mechanism both obser- 
vationally and theoretically, but definitive observations of the smoking-gun are still 
missing. Indeed for all coronal heating mechanisms there are still issues that need 
to be addressed before a definitive conclusion can be reached. For instance, the 
inability to directly 'observe' coronal heating has meant that proxies have had to 
be found. However, it is very difficult using these proxies to distinguish absolutely 
between wave heating or reconnection mechanisms. For example, the drifting of 
heating layers in phase mixing/resonant absorption means that such a wave-based 
mechanism will actually be impulsive in nature and so appear very much like a 



nanoflare heating event ( Ofman et al. 1998 


Moriyasu et al. 


2004 


Klimchuk 


200(5 


Mendoza-Briceno & Erdelyi 


2006 


I. Recent studies by Antol 


in and co-authors ( 


An- 



tolin et al. (2008 2010)) investigate the different observational signatures of coronal 
heating by either nanoflares or the dissipation of Alfven waves. These authors sug- 
gest that coronal rain could be a marker for coronal heating mechanisms, shedding 
light not only on the spatial distribution of heating events in coronal loops but 
potentially also on the actual heating mechanism itself. 

Additionally, on the theoretical front it is clear that simplified models are not 
sufficient for comparison with observations. The response of the plasma in the solar 
atmosphere is very complex and interpreting responses of emission lines is not 
simple either. Forward modelling of theoretical experiments which have physically 
realistic effects included is essential, but the analysis of these results must be done 
with considerable care to determine which are important effects and will take time 
to complete. In the meantime, the argument over which mechanism(s) heat(s) the 
corona will continue. 
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